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Abstract 



We study in more detail the properties of the generalized Beth Uhlen- 
beck formula obtained in a preceding article. This formula leads to a sim- 
ple integral expression of the grand potential of any dilute system, where 
the interaction potential appears only through the matrix elements of the 
second order Ursell operator U2. Our results remain valid for significant 
J> | degree of degeneracy of the gas, but not when Bose Einstein (or BCS) con- 

densation is reached, or even too close to this transition point. We apply 
them to the study of the thermodynamic properties of degenerate quan- 
tum gases: equation of state, magnetic susceptibility, effects of exchange 
between bound states and free particles, etc. We compare our predictions 
to those obtained within other approaches, especially the "pseudo poten- 
tial" approximation, where the real potential is replaced by a potential 
with zero range (Dirac delta function). This comparison is conveniently 
made in terms of a temperature dependent quantity, the "Ursell length" , 
which we define in the text. This length plays a role which is analogous 
to the scattering length for pseudopotentials, but it is temperature de- 
pendent and may include more physical effects than just binary collision 
effects; for instance, for fermions at very low temperatures, it may change 
sign or increase almost exponentially. As an illustration, numerical results 
for quantum hard spheres are given. 
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1 Introduction 



The use of quantum cluster expansions was introduced in 1938 by Kahn and Uh- 
lenbeck Q , who generalized to quantum statistical mechanics the Ursell functions 
U q defined by this author in 1927 |2j. The major virtue of cluster expansions is 
that they provide directly density expansions for systems where the interaction 
potential is not a necessarily small perturbation; in fact, it may even diverge at 
short relative distances (hard cores for instance) while usual perturbations theo- 
ries generate power series in the interaction potential, where each term becomes 
infinite for hard core potentials. Starting from a quantum cluster analysis, the 
Beth Uhlenbeck formula || |], |5|] gives an explicit expression of the first terms of 
a fugacity expansion (or virial expansion) for the grand potential of a quantum 
gas. The expression is valid for any potential, the latter being characterized by 
its phase shifts in a completely general way. 

One should nevertheless keep in mind that the words "density expansion" have 
a double meaning in this context. In a dilute gas, there are actually two dimen- 
sionless parameters which characterize "diluteness" : the product n 1//3 6, where n 
is the number density of the gas and b is a length characterizing the potential 
range (diameter of hard cores for instance), and the product n 1 ^ 3 ^, where At is 
the quantum thermal wavelength of the particles. The former parameter is small 
if, classically, a snapshot of the system shows particles among which almost all are 
moving freely, while the few that interact are engaged in binary collisions only; 
the latter, purely quantum in nature, is sometimes called the quantum degeneracy 
parameter, and remains small provided there is little overlap of the quantum wave 
packets. The validity of the Beth Uhlenbeck formula, as all fugacity expansions, 
therefore requires two independent parameters to be small. 

In a previous article JJ, we discuss one method which conveniently treats 
the two parameters separately, and allows one to include the effects of statistics 
by exact summations while limiting the expansion to the lowest orders in n l ^b. 
The technique is based on the use of Ursell operators U q generalizing the Ursell 
functions (for a system of distinguishable particles), coupled with the exact cal- 
culation of the effect of exchange cycles Ci of arbitrary length I; for short we call 
it the technique of U-C diagrams. The result is another expression of the grand 
potential, which is no longer a fugacity expansion since it includes a summation 
over all sizes of exchange cycles so that statistical effects are included to all orders. 
Truncating the expansion to its first terms (lowest q values and/or low order in 
a given U g ) gives results which remain valid for "dilute degenerate systems" fTj, 
that is for all systems where the potential range is sufficiently small, but where 
the degeneracy of the system may become significant^. In this article we will 

1 For bosons, Bose Einstein condensation is excluded since it requires a summation over an 
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start from an expression of the grand potential which is limited to the first order 
correction in the second Ursell operator [/ 2 , obtained in as a trace over two 
particles of a product of operators. We reduce the trace to an explicit integral 
where the effects of the interactions are contained in a simple matrix element. 
The range of validity of our result is actually similar to that of the calculations 
based on the use of pseudopotentials [§], another approach where the final results 
automatically include the summation of an infinite perturbation series in terms 
of the initial potential. The two methods are comparable, but we think that the 
U-C diagram method provides more general and more precise results, basically 
because it includes the short range correlations between the particles, and be- 
cause all scattering channels with given angular momentum as well as their exact 
energy dependence are included instead of only one constant scattering length^. 
Another point of comparison is the class of methods, for instance discussed in [|10| 
or [IT|, where a renormalization procedure is used in order to obtain expansions 
in terms of the scattering T matrix instead of the interaction potential itself V; in 
the calculations discussed in the present article, no renormalization of this kind 
is needed since, roughly speaking, it is already included in the U q 's, which are 
the building blocks of our method. Nevertheless, as we will see, our method is no 
longer valid when the gaseous system is brought too close to a phase transition 
(superfluid transition of single particles for bosons, of pairs for fermions). 

We begin this article with a study of the expression of the grand potential, and 
show how it can be expressed as an expression that is similar to the well known 
Beth Uhlenbeck formula; actually it can be obtained from it by two simple sub- 
stitutions. We then discuss the physics contained in this general result, as well as 
the changes introduced by the possible occurrence of bound states. In particular 
we consider the effects of exchange between bound and unbound particles, an 
effect which is not contained in the usual Beth Uhlenbeck formalism; this kind of 



exchange may play some role in clouds of laser cooled alkali atoms |L2], [13], nj] for 
which the potential is sufficiently attractive to sustain a large number of bound 
states. In section |3] we apply these results to spin 1/2 particles. Finally, in section 
[|, we discuss the appropriate quantities in terms of which one should describe the 
effects of the potential on the physical properties of the system, and introduce 
for this purpose the so called "Ursell length", which plays a role similar to the 
scattering length a. In the theoretical study of quantum gases, and as already 
mentioned, one frequently used method is to replace the real interaction potential 
between the particles by a "pseudopotential" that has no range (a Dirac delta 
function of the space variables), and to treat this potential to first perturbation 
order; in other words one ignores the distortion of the many-body wave functions 
at short relative distances and the associated effects of the inter particle corre- 



infinite number of interaction terms, a question which we will study in a forthcoming article 

2 This does not mean that one could not improve the theory of pseudopotentials to include 
all phase shifts, since a general expression of the pseudopotential is given by Huang in ||, but 
to our knowledge this has not been done explicitly. 
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lations. The justification of this approach is based on the physical expectation 
that, for a dilute gas, all the effects of the potential should be contained in the 
binary collision phase shifts associated with the potential, which can easily be 
reproduced to first order by a pseudopotential; meanwhile all detailed informa- 
tion on the behavior of the wave function at short relative distances can safely 
be discarded. Our formalism allows one to explicitly distinguish between short 
range effects ("in potential effects") and asymptotic effects (out of the potential), 
which naturally leads to a discussion of this ansatz. The interactions appear in 
terms of a matrix element of an Ursell operator, which depends on the potential 
but does not reduce to it; for fermions at low temperatures, the matrix elements 
contain physical effects which are not included in usual treatments of normal 
Fermi gases. 

2 The grand potential 
2.1 Notation 

The basic object in terms of which most physical quantities will be written in 
this article is the second Ursell operator U 2 , defined by: 

U 2 (l, 2) = [e-^ 1 ' 2 ) - e-fl*i(i)+*i(2)]] (i) 

where -Hi(l) and Hi (2) are single particle hamiltonians, containing the kinetic 
energy of the particle and, if necessary, its coupling to an external potential, and 
where: 

H 2 (l, 2) = H 1 (1) + Hi(2) + V mt (l, 2) (2) 

is the hamiltionan of two particles, including the mutual interaction potential 
V int (l,2). Depending of the context, it may be more convenient to use the sym- 

S A 

metrized operator U 2 ' : 

£^(1,2) = U 2 (l,2) [1+V 2 Pex] = [1+ f ex] U 2 (l,2) (3) 

where P ex is the exchange operator between particles 1 and 2 and rj has the 
value +1 for bosons, —1 for fermions. Moreover, the "interaction representa- 

S A 

tion version" of either U 2 and U 2 ' , obtained by multiplying these operators by 
e PHi(i) e @Hi(2)^ w jjj a | go k e use f u i- we denote them with an additional bar over the 
operator, for instance: 

F 2 (l, 2) = eW^/^e-^ 1 ' 2 ) - 1 (4) 

§ ji 

This operator and its symmetrized version U 2 ' (1,2) act only in the space of 
relative motion of the two particles; they have no action at all on the variables 
of the center of mass. 
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Finally, if the particles are not submitted to an external potential (Hi contains 
only the kinetic energy), it is convenient to introduce the momentum Pq of the 
center of mass of the two particles as well as the hamiltonian of the relative 
motion: 

H rel = — + V int (R) (5) 
m 

[m is the mass of the particles); then, {7 2 (1)2) can be written in the form of a 
product: 

f/ 2 (l,2)=e-^) 2 / 4m x[[/ 2 (l,2)] re/ (6) 

with: 



[U 2 (l,2)} rel =[e-^-e-^] (7) 

2.2 Approximate expression of the thermodynamic po- 
tential 

We now start from relation (46) of which gives the grand potential (multiplied 
by —0) in the form: 

LogZ=[LogZ\ ig + [LogZ\ int (8) 
where [LogZ]^ is the well known value of LogZ for the ideal gas: 

t L °g z \i 9 = -^ Tr { L °S [ l ~ Vze~ PHl ] } (q) 

= 7/rr{Log[l + 77/]} [) 

and where the correction introduced by the interactions is: 

[LogZ) int = z 2 Tr l>2 {u? A (l,2) [1 + V f(l)] [1 + 17/(2)]} (10) 
In these equations, 

z = e^ (11) 

is the fugacity, (3 = 1/ksT the inverse temperature and \x the chemical potential, 
while / is defined as the operator: 



f= i (12) 



ze 
1 — r)ze~ 

Similar results can be found in the work of Lee and Yang, see formulas (II. 8) and 
(11.23) of ref. [TH]. As mentioned in the introduction, equation (ffH) gives the 
correction introduced by the interactions to the lowest order approximation in 
U 2 ; see for a discussion of the higher order corrections. Using the definition 

oiU S 2 A : 

U S 2 A (1,2) = e^V^?7 2 5 ' A (l,2) = 1 + vPex [e^V^e^ 1 ^ _ y 

(13) 
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as well as the relation ze 13111 [1 + r/f] = f, we can rewrite (|T0| ) in the form: 

[LogZ] int = Tr 1|2 {Uf A (l,2)/(l)/(2)} (14) 

which expresses the correction as the average of the operator U 2 ' (1,2) over 
unperturbed distributions functions /'s. 

2.3 Spinless particles and rotational invariance 



For spinless particles, by making the trace in fll4 ) explicit, we obtain: 



[LogZ] mt = -^ J d 3 h j ^ 2 /(k 1 )/(k 2 )aJ A (|k 1 -k 2 |) (15) 
where V is the volume of the system, At the thermal wavelength: 

and where aff A (k) is defined by: 

a SA {k) = < k | U S 2 A | k >= _^ e ^ 2fe2 / m < k | \U 2 S > A ] _ | k > (17) 
Here: 



Xrp \rp 



- rel 



ki 



is the appropriate variable since U 2 ' (1, 2) does not have any action in the space of 
states associated to the center of mass of the two particles; a^ A (k) is a microscopic 
length, independent of V for large systems^, which we will call the "Ursell length" 
- see section |4.1| for a more detailed discussion and a justification of the numerical 
factors that we have introduced. The correction written in (p~5|) is analogous to 
a first order energy correction due to binary interactions, while (k) plays the 
role of some effective interaction (within a numerical factor). 

For instance, if we assume that we treat to first order in perturbation a pseu- 
dopotentialQ of the form: 

3 This is true since Ul' (1,2) has a microscopic range and since the factor V in ( |T7| ) makes 
up for the normalization factor of the plane waves that occur in the matrix element (note that 
all plane wave kets in our formulas are normalized in a finite volume; hence the absence of 
Dirac delta functions of momenta differences in (|l7|). 

4 The most usual procedure is to treat this potential to first order only, since a naive treatment 
of higher orders may introduce inconsistencies. For instance, in three dimensions, it is possible 
to show that all phase shifts, and therefore the collision cross section, of a zero range potential 
such as (|l9|), are exactly zero; on the other hand, they do not vanish to first order (in other 
words, the Born series for potentials containing a delta function is not convergent). For the 
same reason, for a potential such as (|l9|), the Ursell operator f/ 2 vanishes exactly, while it does 
not if the potential is treated to first order. 

For a more elaborate discussion of pseudopotentials going beyond (|l9|) and including waves 
of higher angular momentum, see H. 
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Veff(r) = — 5(r) (19) 
m 

where a is a scattering length (or the diameter of hard cores), we easily obtain: 

< k | U S 2 A | k >~ --ka [1 + v] (20) 

so that ( |T7| ) shows that, in this approximation, the Ursell length becomes inde- 
pendent of the wave number k. For bosons, the correction to the grand potential 
then becomes: 

[LogZ] int ~-2a^- xg 3/2 (z) (21) 

AT 

with the usual notation: 

93/2{z) = J§y J d ' k f^ z ^ ( 22 ) 

while, for fermions, no first order correction is obtained. These results coincide 
with the first order terms of the well-known results of Lee and Yang; see formulas 
(1) and (4) for J = of ref. Nevertheless, in section |], we discuss the 

validity of this first order approximation in the calculation of the matrix elements 
of U 2 and conclude that, for bosons, ( p0|) is a good approximation while it is not 
necessarily the case for fermions. 

At this point, it is convenient to introduce the free spherical waves | j^\ m > 
associated with the relative motion of the particles, with wave functionsf]: 

<^\j ( kL>=ji(kr)Yr(r) (23) 

(with standard notation; ji is a spherical Bessel function, Y™(r) a spherical har- 
monics of the angular variables of r); if | k > is a plane wave normalized in a 
volume V: 

k>=^|EW'[™]*l^™> < 24 > 

l,rn 

If we insert this equality into (fl7|) and take into account the well-known relation: 

£ \yr(k) 2 = ^ (25) 



we readily obtain the result^: 



S,A 



(*)=$>i + l)[l + »7(-l)'] (26) 



i 



5 To normalize these functions (with a Dirac function of k vectors and Kronecker delta's of I 

(0) 



and m), it would be necessary to multiply all the j k / m (r)'s by factors k^2/ir; this operation 
is not necessary here. 

6 We assume rotational invariance, so that the matrix elements of ^2(1, 2) are diagonal in I 
and m. 
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with: 

ni l )(k) = - 



a^(k) = -^e^<j^\m rel \j { l> (27) 



T 

where [^(l, 2)] reZ has been defined in (|7|); rotational invariance ensures that the 
right hand side of (p^ ) is independent of m. These results can be inserted into the 
integral appearing in fll5|) and provide an expression of the correction to the grand 
potential which is a direct generalization of the usual Beth Uhlenbeck formula to 
gases having a significant degree of degeneracy: 

[LogZ]^ = -J^£(2/+l) [l + rj (-iy]x J d% J fhfQcMQtixcflik) 

(28) 



For comparison, we recall the explicit expression of this formula: 
[Log Z]f n f = £(21 + 1) [l + ,(-!)'] / dke-^-kS^k) 

T , J 



(29) 



(we temporarily ignore possible bound states), where is defined in flip]). As 
already discussed in |J, if we replace each of the two /'s in ( |28|) by their low 
density limit e /3(M-^ 2fc2 / 2m ) ; the integration over d 3 K of a Gaussian function in- 
troduces a factor 8 3 / 2 [tt/Xt] 3 and, after some algebra, we recover (^). In other 
words, the following substitutions are necessary to obtain (^) from the usual 
Beth Uhlenbeck formula []: 



8i(k) =>• —ka$(k) 



In equation (|TI1|), the first of these substitutions amounts to adding the terms in 

vf- 



2.4 Bound states 

The operator [t^fgf may be written as the sum of the contributions of bound 
states and of the continuum: 

[C/ 2 ]f;f = £ | $ n )<$„ I e"^" + continuum (31) 

7 There are several equivalent ways to write the Beth Uhlenbeck formula; for instance, an 
integration by parts allows one to replace the product k5i(k) by the derivative d5i(k)Jdk while 
the coefficient 1/At is replaced by 7r/A^. Under these conditions, the second line of (p0|) becomes 
d5i(k)/dk =>- — [k\T) 2 au\k) / it . In other words, the correspondence between ou result and the 
Beth Uhlenbeck formula depends on the way the latter is written. 

The two functions kSi(k) and — k 2 a^ (k) are not necessarily equal but, when multiplied by a 
Gaussian function e~^ h k / m , have the same integral over d 3 k. 
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where the | $ n ) 's are the kets associated with the eigenstate of the relative motion 
of the particles with (negative) energy — E n (with appropriate symmetry for the 
statistics of the particles). It is therefore not difficult make the contribution of 
bound states in fll5|) explicit, which provides the following term: 



[LogZ] 



bound 
int 



(2ir)- b V 2 / d 6 K 



A/(y + q)/(y-q)x 



E i fa i i 2 e 



l3(E n +h 2 q 2 /m) 



(32) 



where ; because these kets have a finite range, and because the plane waves 
| q > are normalized in a macroscopic volume V, the product V | (q | $ n ) | 2 is 
independent of V in the thermodynamic limit, as necessary to obtain an extensive 
correction to the grand potential. If we rewrite the integral of in the form: 



d 3 K / d 3 ge -^ K2 /4m 



E 



we may distinguish between two contributions in the correction: 

Tt rvibound r T rvibound . r T rvibound 



(q | $n) I 2 e^" 
(33) 



(34) 



The first contribution is obtained by ignoring in (|33D the Ty/'s inside the brackets, 
which allows one to integrate over d 3 K; using the closure relations over the plane 
waves | q) and the normalization of the bound states | $ n ) then provides the 
following result for this first contribution of the bound states: 



[LogZ]^ = 2 3 / 2 [A T ]- 3 V^ 



(35) 



which is identical to the term corresponding to bound states in the usual Beth 
Uhlenbeck formula. The second contribution arises from the rest of the product 
of the brackets and is equal to: 



[LogZ] 



bound 
stat. 



(2tt)- 6 V 2 



d 3 K / d 3 ge-^ 2K2 /^ 



/(? + q)/( 



K 



q 



q 



E 

n 



(q I $n) r e 



2 /3E n 



(36) 

It accounts for exchange effects between bound states and continuum states^, 
which explains the appearance of the scalar product (q | $ n ) of a free plane wave 

8 The exchange effects between bound states themsefves will be investigated in another article 
with the study of pair condensation (BCS condensation for fermions); they are higher order in 

Uo. 
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and a bound state wave function; the factors " correspond physically to the 
Boltzmann distribution of the populations of the bound states. 

If we assume that the momenta in the continuum, which have values that are 
limited by the presence of the /'s under the integral, are much smaller than those 
contained in the bound states | $ n )'s, we can replace the product (q | $ n ) by 
(q = | $ n ); this shows that the effect of statistics is more important for bound 
states with a wave function with a significant integral over space, as for instance 
the ground state wave function which has no node; states with rapidly oscillating 
wave functions give almost no correction to the usual Beth Uhlenbeck formula. 
Generally speaking, for bosons, the effect of statistics is always to increase the 
Beth Uhlenbeck term (]35|) . For fermions, the situation is more complicated: if the 
gas is only slightly degenerate one has f 2 < f so that the terms in 77 dominate in 
(|36|), leading to a decrease of the effects of bound states; but if the gas is strongly 
degenerate, there seems to be no general rule, and exchange of bound states with 
the continuum may either enhance or reduce their contribution. 

2.5 Spins 

An easy generalization is to include spins; the only difference is that all the traces 
must now also include spin states. It is then convenient to replace the functions 
/(k) of momentum by spin operators /s(k) (corresponding to two by two matrices 
for spin 1/2 particles) which are functions of k and are defined as: 

ze' 13111 

(m s I / s (k) I m s ) = (m s , k | - _ _ pHl I m s , k) (37) 

where Hi is the one particle hamiltonian, including kinetic energy as well as 
coupling of the spins to the magnetic field (if the particles carry magnetic mo- 
ments). Similarly, because now the exchange of particles must also include their 
spin states, a^ A (k) becomes an operator Es(A;) which acts in the space of the 
states of two spins: 

00 

= X> + 1 ) <$(*) [1 + V(-1) 1 PI] (38) 
1=0 

where Pf x is the exchange operators of two spins. This leads to the following 
generalization of fll5l) : 

[LogZ]^ = J d% J d 3 k 2 Tr 5li52 {fs(ki)f s (k 2 ) S 5 (A;)} (39) 

3 Spin 1/2 particles 

We now apply the preceding calculation to the study of the magnetic suscepti- 
bility of a dilute gas of fermions or bosons with spin 1/2 (spin polarized atomic 
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hydrogen provides an example of spin 1/2 bosons |17|, [lBfl). We assume that the 
one particle hamiltonian is: 



-a. 



(40) 



2m 2 

where Uq/2tt is the Larmor frequency in the (homogeneous) magnetic field and 
a z the Pauli matrix (operator) associated with the component of the spin along 
the field; we set: 

which gives the following values for the matrix elements[]of /s(k) defined in 



/±(k) = 

For short we will write: 

/±(1) = /±(ki) 
With this notation we have: 
V 



z±e 



-/3h 2 k 2 /2m 



1 — 7]z±e 



-j3h 2 k 2 /2m 



/±(2) = /±(k 2 ) 



(42) 



(43) 



[LogZ] 



'.</ 



V 



8tt 3 



d ki Log 



V 



1 — 7]z + e 



-f3h 2 k 2 /2m 



1 — r]z_e 



-f3h 2 k 2 /2m 



d 3 hLog [l + 77/ + (l)] [l+77/_(l)] 



(44) 

For a gas at equilibrium the operators /s(k) are diagonal in the basis correspond- 
ing to a quantization axis parallel to the magnetic field; we then have: 



TWJ/s(ki)/ s (k 2 ; 



as well as: 



Tr Sl {/ s (k 1 )}xTr 52 {/ 5 (k 2 )} 
[/+(!) + /-(!)] [/+(2) + /-(2)] 



Tr 5li52 {/ 5 (k 1 )/ s (k 2 )PfJ = /+(l)/+(2) + /_(l)/_(2) 



(45) 



(46) 



(the latter result arises because the trace gets non zero contributions only from 
the two spin states | +, +) and | — ,— ), which are invariant under the effect of 
Pex)- We therefore have: 



[LogZ] 



hit 



(2vr) c 



d'h I d A k 2 la^ik) [/ + (l)/ + (2) + /_(l)/_(2)] 



+2a v (k) /+(!)/_ (2) 



(47) 



9 We remind the reader that, except for an ideal gas, the /±'s are not the populations of the 
one body density operators, but differ from them by density corrections [E5[. 
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where: 

k=^M (48) 

while afj A (k) is defined by ( |TTD - or equivalently ( p6| ) - while au(k) is the un- 
symmetrized version^ of the Ursell length: 

M*0 = < k I C/ 2 I k >= V(2/ + l)a^(A;) (49) 

(it would correspond to distinguishable particles). If, as in the beginning of 
section we treat the pseudopotential (M) to first order, for fermions we 
obtain a^{k) = and ajj{k) = a; the first order correction is now given by: 

[LogZ] int ~-2^ xg 3/2 (-z) (50) 



which coincides with the result of Lee and Yang (equation (1) of [16| for J = 1/2); 
but, again, we note that a critical discussion of this first order calculation is made 
in section |] which shows that, for fermions, the results of the Ursell approach 
may be different from those of a pseudopotential theory. 



3.1 Density of particles and of energy; magnetization. 

The number density of the gas is obtained from the relation: 



n 



V dz 
Similarly, then density of internal energy is given by: 



V- 1 z— [LogZ] (51) 



^ = ^ = -V" 1 ^[LogZ] (52) 



while the "magnetization"!^]] is equal to: 



M = hih u ' iZ] (53) 



10 Note the factor 1/2 which does not appear in (|Tj); we choose this convention since the 
same factor 1/2 appears in the definition ([|) of the symmetrized version of U2', in this way, if 
exchange effects are ignored (high temperature limit for instance), the various Ursell lengths 
become equal. 

n What we call here magnetization is not a real magnetic moment (ampere square meter) 
but a dimcnsionlcss number equal to the sum of the average values of a z of all atoms; in other 
words, the maximum value of M (complete spin polarization, all spins parallel) is equal to the 
total number of particles. 
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We therefore have to vary either z, or (3, or uiq in formulas ( f44| ) and (f47|); we can 
then use the simple relations: 



and: 



as well as: 



dz± = z± ( — ± —^-d(3 ± Ph—^- 



df± 
dp 



dz-i 



U [i + vf±] 



h 2 k 2 hjo 



2m 



± 



(54) 

(55) 
(56) 



We then obtain, by making use of the symmetry in the indices 1 and 2: 

n =(2^1 d 3 Aa [/+(!) +/-(!)] + 

J d 3 h J d 3 k 2 ^a^ik) U(l)f + (2) [1 + 77/ + (l)] +id. (f + & f. 

+2a v (k)f + (l)f-(2) [2 + V f + (1) + tj/_(2)] 



(2^ 



(57) 

where id. (/+ -v=> /_) symbolizes the same expression where /+ and /_ are inter- 
changed. The calculation of M is almost the same, except that now, because of 
(|54"D , /+ and /_ introduce different signs; the result is: 



^ =(2vr)- 3 [/+(!) -/-(!)] + 

A j 1 



(27T) 



| d 3 h J d 3 k 2 |2ag A (A;)/ + (l)/ + (2) [1 + V f + (l)} -id. (/+ ^ /_) + 

+2i70tr(*)/ + (l)/-(2) [/ + (1)-/_(2)]| 

(58) 

Finally, the calculation of the internal energy w provides the result: 

w = w E + w M (59) 

where we is the density of energy associated with the external variables of the 
particles (kinetic and potential energy): 



w E 



:(2tt) 

A ji 

(2^y 



ft 2 *; 2 



d 3 ^ [/+(!) + /_(!)] 



2m 



J d% J d 3 k 2 |2a5 A (fc)/ + (l)/ + (2) 

.... . , . , , r r , ^ 2 fc? 



2m 



+ id. (/+ ^ /_ 



+2au(k)f + (l)f42) 



h 2 k 2 2 

2m 
(60) 



[l + r?/ + (l)]^ + [l+7 7 /_(2)] 
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while wm is the density of magnetic energy: 



w M 



ftw M 



(61) 



We now assume that the magnetic field is sufficiently low so that the preceding 
expressions may be limited to their first order expansion in uq ~ 0. We then use 
the relations: 



, , . df± dz± f ,9f (3ftw 

f±-f + ^— ^—^o = / ± z- — = / 

oz± oojq oz 2 



i ± ph^- [i + nf\ 



(62) 



which provide: 

n{u ) ~ n(0) = (47T 3 )- 1 / d%f(l) 



(2^y 

,S,A, 



d 3 h / d 3 k 2 f(l)f(2)x 



x 4 [1+77/(1)] [a^(k) + au (k) 
and for the zero field magnetic susceptibility of the gas: 



(63) 



X 



M 
du V 



(3h / d 3 hf(l) [1 + 77/(1)] 



(27T) 



/ d 3 fc! / d 3 fc 2 /(l)/(2)x 



x <{ 2 a5 A (fc) [1 + 7^/(1)] [2 + 27^/(1) + 77/(2)] + 



+ au (k) [[/(!)- /(2)] 2 + 277/(1) [1 + 77/(1)]]} 



(64) 



3.2 Physical discussion 

Equations fl5T|) and (|58| ) show how quantum statistical effects modify the equation 
of state and the magnetic susceptibility of a gas. What determines the properties 
of the dilute system is not directly the matrix elements of the potential Vi nt , but 
those of the operator {7 2 - Of course, if the potential is weak, the calculation 
may be limited to first order perturbation theory, so that U2 may be replaced 
by — PVint] the two operators then become equivalent. But, in more realistic 
situations, for instance in atomic systems where the interaction potential at short 
relative distances becomes very large, the effects of the potential cannot be treated 
to first order. An obvious difference between C/ 2 and Vi nt is that the former 
depends on one more characteristic length than the potential, namely the thermal 
wavelength At, so that one can expect that its range will depend explicitly on 
the temperature. It is sometimes argued that a description of the interactions 
in terms of the scattering length is sufficient at very low temperatures, because 
the range of the potential Vi nt remains always much smaller than the de Broglie 
wavelengths of the particles. As far as the matrix elements of U2 are concerned, 



14 



this is of course true if the interaction potential is treated to first order only but, 
precisely at very low temperatures, the higher order terms become significant and 
TJi may acquire a range that exceeds by far that of the potential itself. We will 
see examples of this phenomenon in section |j. 

We now discuss the limit of validity of the results obtained in the preceding 



section or, since they are equivalent, of equations (|T4]) and (0). If the particles are 
bosons, our theory is limited to gases which are not too close to Bose Einstein 
condensation - while of course the degeneracy may be much more pronounced 
than for the Beth Uhlenbeck formula since z does not have to be small. The reason 
is that, when a system of bosons approaches the region of quantum condensation, 
the distribution function / starts to build up a singularity at low energies. One 
can then see that, when z — > 1, the partition function becomes more and more 
sensitive to terms of higher and higher order in U2, in U3, l7 4j etc. , since they 
contain larger and larger number of functions / (or of factors 1 + rjf). In other 
words, terms which normally remain small corrections become dominant when 
the point of Bose Einstein condensation is approached; see the discussion already 
given in ]7j , where it is emphasized^ that a theory limited to first order in {7 2 
would predict the disappearance of the Bose Einstein condensation phenomenon 
and its replacement by a simple crossover between two regimes. 

For fermions, the discussion is different since no special phenomenon takes 
place when z reaches one; when the gas is cooled at constant density, this merely 
corresponds to a cross over region where the gas is becoming degenerate. When 
the temperature is decreased even more, a stronger degeneracy builds up while 
z =er^ becomes larger and larger (/3 increases while \x remains almost constant 
if the density is fixed). It is therefore clear that z itself can not be an expansion 
parameter in this region; but our perturbation series is not a z expansion and it 
may still converge for degenerate systems, provided the matrix elements of U2 
are sufficiently small. In the next section, we will see that these matrix elements 
are equal to the product A^a^, where At is the thermal wave length and a v is 
some microscopic length (the Ursell length) that we will define more precisely 
later. Let us for instance discuss fl5"T|) where the interaction corrections involve 
a double integral over two momenta. The first introduces the number density n; 
the second contains a function /(l — /), which in the limit of low temperatures is 
non zero only in an energy slice of width (3" 1 at the surface of the Fermi sphere; 
the corresponding width Ak in terms of momentum is given by: 

= r 1 (65) 

m 

(where fikp is the momentum at the Fermi surface) so that the result of these 



12 Unfortunately, in j?J an assumption of this discussion is not made explicit, namely the 
fact that the matrix elements of XJi should be positive (dominant character of the attractive 
interactions in the matrix element). 
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integrations is: 

n (A^ajj) x Ankp- — -j- = 2nk F ajj (66) 

The small parameter of the expansionP] is therefore the product kpciu, which is 
not very different from n l ^au for degenerate gases, and indeed remains small if 
the average distance between the particles is much larger than ajj. This result 
is not surprising: indeed, in the usual theory of dilute non ideal Fermi gases 
pTOfl , the expansion parameter is independent of the temperature, which plays no 



particular role as long as the system is strongly degenerate. 

One might be tempted at this point to conclude that the generalized Beth 
Uhlenbeck formula is valid for arbitrary degeneracy of a fermionic system; never- 
theless, in the next section, we will see that the factor ajj itself may become very 
large at very low temperatures, which in turn increases the expansion parameter 
and, automatically, limits the validity of the expansion at some point. We will 
assume that this somewhat unexpected phenomenon is a precursor of BCS pair 
condensation; if this is the case, the low temperature limit of the validity of the 
generalized Beth Uhlenbeck is that the gas should remain a normal Fermi gas 
and, even, not to be too close to condensation. In general, for bosons as well as 
for fermions, we can therefore conclude that our calculations remain valid as long 
as Bose Einstein or BCS condensation is not too close. 



4 Matrix elements of JJi 

We now study in more detail the values of the basic ingredient that we use 
to describe the effects of the interactions on the thermodynamic properties of a 
quantum gas, namely the diagonal matrix element of U 2 ' - together with those of 
the un-symmetrized operator if the particles have spins, see section |3|. A natural 
question then is the following: to determine these coefficients, is it sufficient 
to characterize the potential in terms of its binary collision phase shifts, which 
determine the asymptotic behavior of interacting wave functions, or is it also 
necessary to include some information on the behavior of the wave functions at 
short relative distances, inside the potential? For bosons at very low temperatures, 
is it possible, even more simply, to reason in terms of the scattering length only? 

We have already mentioned in the introduction that, in the literature on low 
temperature dilute systems, it is often considered as physically obvious that the 
short range properties of the potential play no role in the thermodynamics; it is 
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We note in passing that the phase occupation factor (1 — /) plays an essential role in this 



argument. This factor occurs in all expressions, for instance, in (58), the bracket [/+(1) — /-(2)] 
may be written as [/+(1) — 1 + 1 — /_(2)]. If we had ignored this factor, we would have found 
nau^r as the expansion parameter, which is a temperature dependent factor (as in the usual 
Beth Uhlenbeck formula). 
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then possible to replace the real potential by an effective potential (or pseudopo- 
tential) which has the expression already given in (|19|) in terms of the scattering 
length a. With a potential having zero range, it becomes of course meaningless 
to take into account the distortions of the wave functions inside the potential, 
in other words the contribution of particles "in the middle of a collision". But 
with real potentials, it is not obvious that these contributions are non existent, 
in particular since statistical effects arising from particle exchange are expected 
to be more important when the particles are close. We discuss here this question 
in the context of the generalized Beth Uhlenbeck formula that we have obtained, 
with the help of simple examples such as a step like potential mimicking a real in- 
teraction potential, or even hard spheres. We first discuss the Ursell length au(k), 
which we have already introduced above, and which depends on the temperature 
while a does not. 

4.1 Ursell length 

If the interaction potential V int were sufficiently weak, we could use a first order 
perturbation expression of the Ursell operator: 

U 2 = - /Ve-^ VW^'-^ + ... (67) 
Jo 

which would lead to the following expression: 

< k | U 2 | k >= -p < k | V int | k> +... (68) 
Inserting (|i~9|) into this result provides: 

< k U 2 k >= -P^rr- = ^ 69 

where V is the macroscopic volume and the thermal wavelength flTEp. 

By analogy with this first order calculation, it is convenient to characterize 
the diagonal matrix elements in terms of a length, au{k), which we call the Ursell 
length and which we have defined in fl49|) as: 

a u{k) = -X < k | U 2 | k > (70) 

2 A rp 

This very definition ensures that, within the theory of the pseudopotential, au(k) 
is exactly equal to the scattering length a. This length therefore provides a conve- 
nient tool for a discussion of the validity of this approximation: as long as ajj(k) 
remains very close to a for all relevant values of k, the theory of pseudopoten- 
tials and ours provide strictly equivalent results. For identical particles, what is 
needed is the symmetrized Ursell length, already defined in ( |TTD and Q26]), which 
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allows us to write the correction to the grand potential for spinless particles in the 
form fll5|). This equation ([TJ) expresses that the symmetrized Ursell length gives, 
within the numerical factor in front of this integral, the crossed contribution of 
two velocity classes ki and k 2 to the grand potential of the system (its pressure). 
If the particles have spins, we have to use (|39| ) and (|38|). 

In section | 2.3| , we introduced spherical orbital variables from the free spherical 
>. At this stage, it is convenient to use spherical waves that are 
normalized in a sphere of volume V = 47ri? 3 /3. We therefore introduce a new 



waves | f' 



notation, | <fifi m >, but these kets are simply proportional to the 

I ,/,(°) ^- rr I M ^ 
I Yklm x kl | Jklm ^ 



■ (0) 
■3 klm 



>: 



(71) 



where x^ is a normalization coefficient which in the limit of large values of the 
product kR is equal to: 



Xkl 



2k 2 



Introducing in the same way the normalized kets 
interaction potential, we can write: 



(72) 

fk'im > i n the presence of the 



2vr 
X 



2 l x kl\ 



k' 



(73) 



where e(k) and e(k') are the energies of the free and interacting states. 

In what follows we discuss the values of the Ursell length and its k dependence 
with the help of a few examples. 



4.2 Steplike potential 

In order to simplify our discussion, we now consider an interaction potential 
made of a hard core of diameter xb (with x < 1) and of an attractive part 
(Vint — —Vq) from relative distance r = xb to r = b; see figure 1. Our discussion 
is in fact more general, but this kind of simplified potential is a convenient way 
to mimic the effects of a more realistic interaction potential, containing strong 
repulsion at short distances and Van der Waals attraction at large distances. For 
instance, it is not difficult to find a relation between Vo, x and b which ensures 
that the scattering length of this potential vanishes (compensation of the effects 
of attraction and repulsion at low energies). Does this imply that the matrix 
elements of U2, that is the Ursell length ajj{k), also vanish? Not in general, since 
the scalar products < (pf] m \ (p k > l m > in ( [73]) are not only sensitive to the changes 
of the interacting wave functions outside of the potential (which do not occur if 
the scattering length vanishes) but also to their values inside the potential. In 
other words, for a degenerate dilute gas, all the effects of the potential are not 
necessarily contained in the scattering phase shifts. 
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We nevertheless note that, when the range b of the potential tends to zero, 
the corrections to the scalar products are necessarily of third order in b (at least), 
while first order terms occur in the phase shifts (for instance, at very low temper- 
atures, the Beth Uhlenbeck formula predicts a correction which is proportional 
to 6Ay). Therefore, if the potential range is sufficiently small (compared to the 
two other microscopic distances, the average distance between particles n -1 / 3 and 
Xt), one may limit the calculation to first order in b so that the contribution of 
the distortion of the wave functions inside the potential may be ignored. Within 
this approximation, all effects of the potential on the thermodynamic properties 
are indeed contained in the collision phase shifts. A brief similar discussion in 
terms of "shape dependent terms" (as opposed to phase shift dependent terms) 
appears in the two last paragraphs of section 5 of ref. ||T9| 

Consequently, in all calculations of the thermodynamic quantities which are 
limited to first order in the potential range, it is sufficient to characterize the 
potential by its long distance effects on the wave function (phase shifts) only, while 
short range effects are irrelevant; if all collisions take place at very low energies, 
all effects are then contained in the scattering length only. On the other hand, 
this is not necessarily the case if higher order terms in the potential range are 
included, and corrections which originate "inside the potential", corresponding 
to the contribution of particles "in the middle of a collision", may appear; see 
also [[HJ. Actually, the question remains open as to whether is would be possible, 
by some mathematical transformation, to express these "in potential effects" in 
terms of the phase shifts only; this is for instance possible in the absence of the 
phase occupation factors, since one then gets the usual Beth Uhlenbeck formula 
for which such a transformation is known. In the presence of these factors, we 
have made efforts to investigate the possibility that a similar simplification takes 
place, but we have not been able to prove it. A possible conjecture is that short 
range and long range effects are in general independent from each other, but this 
remains to be proved by a precise example. 



4.3 Hard spheres 

For hard spheres, the interacting wave function does not penetrate into the poten- 
tial; this is a special case where all physical effects of the potential are necessarily 
contained in the collision phase shifts. The potential range b coincides in this case 
with the scattering length a. But this is not sufficient to ensure that the theory 
of pseudopotentials should be equivalent to our results, and we now discuss this 
question. 



4.3.1 Analytical calculation 

In order to obtain the correction to the thermodynamic potential for hard spheres, 
we now perform a calculation of the scalar products which appear in ([731). 
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s waves For free waves, the quantification for the wave numbers is given by: 

k = n i < 74 > 

while, for hard spheres of diameter a, the interacting waves satisfy the relation: 

/.' = ^— (75) 

then become 



n tt 
R-a 



where n and n are integer numbers. The scalar products of 
functions of these numbers: 



< <4°o,o I (pk',o,o >z 



R 



dr sin 



mrr\ 



R 



sm 



n'n (r — a) 



R-a 



(76) 



which, after a simple calculation, leads to: 



< v4°o,o I ft',o,o > 



rnra 



it sm 



An 2 [1 - 1] 



R 



n 



n (1 - a/R)f [n + n (1 - a/R)} 2 



(77) 

This result, when inserted into (|73|) for Z = 0, provides the contributions of s 
waves to the Ursell length. Equivalent results for infinite volume (only s wave) 
are given by Lee and Yang fl2T| ; we have checked that our numerical results are 
compatible with those of these authors. Since na/R = ka/n, equation ([771) shows 
that the scalar product is peaked around a value of n given by: 



n 



n - E(-) 



(n 



where E is the integer value; for low energies, b C 1 and the preceding equation 
reduces to n = n. 

Assume for instance that we are interested in calculating the Ursell length to 
first order in a, and in the s wave channel only. Equation (|77|) shows that only 
the term n = n contributes and that the corresponding value is: 



mxa 



sm 



4n 2 [1 - |] 



so that 



R p*i]V[2-£]- 

becomes to this order: 



(79) 



aP(k) 



2?r R 
A| 2P 



n 2 h 2 n 2 



rn 



{R-a) 



1 

R 2 



(80) 



The first order value is therefore merely equal to a; but we will see that a^\k) 
may strongly differ from this value if higher orders effects in the potential range 
are included. 
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larger / values Formula ([76]) can be generalized to non-zero angular momen- 
tum /. The sine functions of the free and of the hard-sphere wave-functions have 
to be replaced by a spherical Bessel function and a linear combination of this 
function and a Neumann function of order I respectively. After normalization we 
obtain for Irl > a: 



2a 2 k* 



ji(k a)ni(k r) - m(k a)ji(k r) 



R[{ka) 2 m 2 {ka) -a/R] 

(81) 

mi{k a) stands for "modulus" pi], 9.2.17]: 

mi{k'a) = ^jf(k'a)+n 2 (k'a) (82) 
The quantification conditions have to be modified; ( [75] ) now becomes implicit: 



k R + T]i(k ) = nix + 



7T 



i3) 



where rji (k ) stands for the /-wave phase shift for wave number k\ which is given 
by: 

/ ji(k a) 

fji{k ) = arctan — — m k m (84) 



ni(k'a) 



In this definition we limit the values of the inverse tangent to the interval [— | ; | [ 
and m k ' counts the number of times that — by increasing k'a from zero to the 
final value — the value of inverse tangent jumps from ~ to — ~. In this way we 
obtain rji as a continuous function of k' . For large values of k'a, the effects of the 
centrifugal barrier become negligible and the quantification condition becomes 
independent of I so that it reduces to the s-wave expression (]75|). 

The normalized free wave function is given by (0), so that we arrive at the 
following expression for the square of the scalar product: 



(o) I 



2a(ka)(k'a) 
R[{k'a) 2 - {ka) 2 }\ {k'a) 2 m 2 (k'a) - a/R 



15) 



4.3.2 Numerical results 



Inserting ([77|) and Q55 ) into ([73]) and then into fl4"9| ) - or ( PB| - provides the Ursell 
length, which in turn determines the diagonal matrix element of U2 - or U 2 ' for 
identical particles. For brevity, we just give the results of our numerical calcula- 
tions of au(k)\ in other words we only discuss interactions between particles in 
different spin states if the particles are indistinguishable. But there is no diffi- 
culty in treating the general case, since the calculation of af/ A (k) is very similar; 
moreover, as soon as the wave number k is sufficiently small, one simply has 
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afj{k) ~ 2ajj(k) and a§(k) ~ 0. The results concerning ajj(k) are shown in fig- 
ures 2 and 3 (more details about the calculations can be found in p2]). There are 
two parameters in the problem, a and At • As convenient dimensionless variables 
we choose, either the product kXr, or ka ; the ratio a/Xr is then kept as a fixed 
parameter, which has small values either when the potential range is very small 
or when the temperature is low. Figure 2 shows that, when a -C At, and as long 
as /cAt remains smaller than 1, the Ursell length is indeed equal to a (with a good 
accuracy), so that the theory of the pseudopotential is fully justified. For higher 
temperatures, we note that the low energy values of au{k) become noticeably 
different from a. This is not surprising since, when T increases, we progressively 
reach a classical regime of small wave packets; they no longer have a much larger 
spatial extent than the potential, while this is necessary for the approximation 
of the pseudopotential to apply. We therefore concentrate on low values of the 
ratio a/ At- 

If the system is made of bosons, since the distributions /'s contained in (|15|) 
have at most the same width as Gaussian thermal exponentials, the values of /cAt 
that are relevant in the integral are comparable to 1, or smaller; figure 2 then 
shows that the use of our theory or of the theory of pseudopotentials leads to the 
same results; see the discussion of section [2~3| , in particular formula (pl|). But 



assume now that the system is made of fermionsf^], and that the temperature is 
progressively lowered at constant density. If the system is degenerate, the width 
of the functions /'s is determined by the Fermi momentum kp, which in turn 
depends on the density; in other words the width is practically independent of 
the temperature. Meanwhile, if the temperature is lowered more and more, At 
increases so that the product k\j< can take arbitrarily large values inside the 
integral. On the other hand, ka remains limited to values smaller than 1 since, 
for degenerate gases where the distance between the particles is comparable to the 
inverse Fermi momentum in a degenerate gas, the diluteness condition n 1//3 a <C 1 
yields: 

1 > k F a > ka (86) 

The departure from 1 of the curves of figure 2, which fall well below this value 
when /cAt increases, shows that significant discrepancies from a pseudopotential 
theory may indeed be obtained. This is even more visible in figure 3, which uses 
a different variable, the product ka: even for small values of the ratio o/At, sig- 
nificant departures of the Ursell length from a are obtained for small values of 
ka; in other words, fl50|) is no longer a good approximation of the generalized 
Beth Uhlenbeck formula. Actually, the smaller a/Xx, the more pronounced these 
departures are; we have an illustration of the consequences of the temperature 
dependence of U2, where the effects of the potential are more and more delocal- 
ized by thermal effects while At increases; consequently, variations of the matrix 
element occur even if the range of the potential is very small and even if ka <C 1. 



We are dealing here with the description of interactions of fermions in opposite spin states. 
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In figure 4, we plot the variations of the diagonal element < r | [L^L^ | r > as a 
function of the relative position r; the results show clearly how the second Ursell 
operator acquires a longer and longer range at decreasing temperatures. 

A striking feature of figure 3 is the change of sign of the Ursell length which 
takes place when ka increases from zero. The origin of this change is under- 
standable from (|73|), from which one can convince oneself that the contribution 
of low values of k becomes dominant as soon as k is sufficiently large; this is 

because, while the scalar product < <ffi m | W \ m > with k' ~ decreases rel- 
atively slowly when k increases (as k~ 2 ), the exponential e / 3e ( fc )-/ 3e ( fc ) varies much 
more rapidly and so that it makes small values of k dominate the sum. In other 
words, what determines the diagonal matrix element of U 2 is the contribution 
of interacting states that have a very small relative energy; because these states 
evolve more slowly in time than the free wave packets of energy h 2 k 2 /m, the net 
differential result is equivalent to an attraction. We therefore come to the conclu- 
sion that even hard cores can result in an effective attraction between fermions of 
opposite spins at the surface of a Fermi sphere, provided that the temperature is 
sufficiently low. But other interesting features also appear; for instance there is a 
strong dependence of the matrix element on the relative momentum and, for some 
value of k, the Ursell length (the effective interaction) vanishes; probably more 
important is the fact that the effective interaction increases almost exponentially 
as a function of k 2 , which is nothing but the square of the relative momentum of 
the interacting particles. This shows that the correction to the partition function 
is dominated by processes that take place preferably between particles having 
almost opposite momenta on the surface of the Fermi sphere (assuming that they 
have opposite spins). 



4.3.3 Validity of the Ursell expansion 

We now come back to the discussion made at the end of section |3.2| concerning 
the validity of the Ursell expansion for a dilute gas of fermions. The result of 
this discussion was that the expansion parameter is the product n 1//3 a;y. This 
parameter would remain a constant as a function of temperature for a gas of 
constant density if the Ursell length did not vary too much as a function of 
temperature, so that the situation would be simple. But in fact, we have actually 
found for Fermi gases that the maximum value of a\j becomes larger and larger 
at low temperatures. This automatically limits the range of validity of the first 
order U 2 theory to temperatures at which the relevant values of the Ursell length 
are not too large. 

To determine a limit temperature, we will look for an asymptotic expression 
of the Ursell length valid for wave vectors of the magnitude of the Fermi wave 
vector k ~ kp. In this case, the terms corresponding to low values of k! in the 
sum over the states of interacting pairs of ( |T3"D are dominant. When k' <^k F the 
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square of the scalar product (|85| ) reduces to: 



(0) I 
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sin kpa (87) 



We replace the square of the scalar product in ([73]) by this expression and trans- 
form the sum into an integral: 



du{kp) oc 



R 2 
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oc - A f sin 2 fc F a e^ (89) 
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where ep is the Fermi energy (proportional to n 2 ^ 3 for a strongly degenerate gas). 
We are interested in the limit n 1 / 3 a <C lor /c]/ 3 a <C 1; we can thus replace the 
sine by its first order approximation, which yields: 

The validity criterion of the Ursell expansion then becomes: 

n 1/3 a—^—r^e^ k - T « 1 (91) 
{kfXt) at 

which shows that the theory is valid only when the temperature satisfies the 
approximate condition: 

kBTc ~ t r xTfx i ^21 ( 92 ) 

Log [nX T (X T /ay\ 

Inside the logarithmic function in the denominator, both the factor n\ T and 
(Xt/ci) 2 are larger than one, especially the latter for a very dilute gas; but since 
the logarithm has slow variations, whatever parameters we choose, in practice 
the limit temperature fixed by this condition is not lower than is about a tenth 
of Fermi temperature, a rather a high temperature compared to the transition 
temperature for Cooper pairing. 

Our conclusion is that, for any type of interaction, and even if the density 
of the gas is fixed at a very low value, at sufficiently low temperature the Ursell 
length becomes larger and larger so that the Ursell expansion is no longer a 
good expansion. The maximum value of the Ursell length is obtained for pairs of 
fermions with opposite spins and maximum relative momenta, which corresponds 
to two fermions having opposite momenta on the surface of the Fermi sphere; this 
is reminiscent of the BCS pairing phenomenon, while the phenomenon takes place 
at much higher temperatures (the right hand side of (|92D does not coincide with 
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the standard expression of the critical temperature) so that stricto sensu it can 
not be called a precursor of this transition. Whether of not it is related to this 
transition, the change of sign of the Ursell length will have a strong effect on 
the two body correlation function in the system (we have seen in p5| that the 
second Ursell operator is closely related to the short range properties of the two 
body density operator). These results are also reminiscent of the well known work 
of Luttinger and Kohn |23[] who predict the occurrence of superconductivity in 
purely repulsive systems, as well of the more recent work of Kagan and Chubukov 
|4}| who also predict in this case p wave superfluidity in a dilute Fermi gas. We 



are planning to investigate this connection in more detail in a future article. 

To summarize, for both fermionic of bosonic systems, the validity of the Ursell 
expansion is limited, even for a very dilute gas, to temperatures sufficiently above 
any transition point. Mathematically, in the case of bosons the divergence of the 
series arises from the distribution function (or operator) /, while in the case of 
fermions its origin is the increasing value of the matrix elements of JJi themselves. 



5 Conclusion 

Our formalism provides a systematic treatment of the interactions in a dilute gas 
where the basic objects are not the matrix elements of the potential itself but 
those of temperature dependent operators. In |p5| |, we have shown how micro- 



scopic, short range correlations between particles could explicitly be taken into 
account and calculated. In the present article, we investigate the macroscopic 
properties of the gas by basing the discussion on expression (fH]), which resem- 
bles a first order perturbative expression of an energy correction, while it actually 
is rather different. This is mostly because the matrix element which appears in the 
expression is not the matrix element of the potential itself, or of some variety of 
pseudopotential, but the matrix element of the second Ursell operator U2, which 
corresponds physically to a local Boltzmannian equilibrium. Except of course in 
trivial cases where the interaction potential is indeed weak for all values of the 
relative distances of the particles, which allows for a first order treatment of the 
potential, this introduces significant differences; the major reason is the temper- 
ature dependence of the matrix elements of U2, which is in general more complex 
than being merely proportional to (3 (as would be the case in usual perturbation 
theory). This matrix element is conveniently described in terms of the Ursell 
length. We have seen in particular that, for fermions at very low temperatures, 
effective attractions at the surface of a Fermi sphere may take place, indepen- 
dently of the repulsive or attractive character of the potential itself. Moreover, 
our formalism contains naturally effects such as the statistical exchange between 
bound molecules and free particles, which may play some role in experiments 
with alkali atoms at very low temperatures [12], O, ff4l. 
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Note added 



After completion of this manuscript, we received a preprint from H. Stein et 



al.[2"6]| who, in the case of bosons, find that the many body scattering length 
undergoes a divergence when the superfluid transition is approached, an effect 
similar to what we obtain for fermions. The significance of this similarity has not 
been examined in detail yet. 



Figure captions 



Figure 1: Steplike potential used to discuss compensation effects between the 
attractive and the repulsive part on the scattering length; Vq is the depth of the 
attractive potential, b its range; a hard core of range xb is assumed (x < 1). While 
a relation between these three parameters can be found to make the scattering 
length a vanish, the effect of the potential on the wave function inside its range 
remains significant. 

Figure 2: Variations of the ratio between the Ursell length au{k) and the 
range a of the hard core potential as a function of the dimensionless variable 
k\x- The ratio a/X T , with values shown in the upper right, is a parameter which 
takes small values if either the temperature is sufficiently low or the radius of the 
hard cores sufficiently small. 

Figure 3: Variations of the same quantity as in Figure 2, but as a function of 
the dimensionless variable ka\ this representation is convenient for the discussion 
of fermions (where it characterizes interactions between particles in opposite spin 
states). 

Figure 4: Variations of the diagonal matrix element < r | [t^],.^ | r > as a 
function of r, for various values of the parameter a/A^; at low temperatures, the 
thermal increase of the range of the operator is clearly visible. 
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Figure 1: Potential used for the discussion in section 4.2 



